arXiv:1506.04466vl [astro-ph.GA] 15Jun2015 


ApJ Letters, in press 

Preprint typeset using style emulateapj v. 5/2/11 


PRE- AND POST-BURST RADIO OBSERVATIONS OE THE 
CLASS 0 PROTOSTAR HOPS 383 IN ORION 

Roberto Galvan-Madrid\ Luis F. Rodri'guez^-^, Hauyu Baobab Liu^ 

Grainne Costigan"*, Aina Palau\ Luis A. Zapata\ Laurent Loinard^ 

ApJ Letters, in press 

ABSTRACT 

There is increasing evidence that episodic accretion is a common phenomenon in Young Stellar 
Objects (YSOs). Recently, the source HOPS 383 in Orion was reported to have a x35 mid-infrared - 
and bolometric - luminosity increase between 2004 and 2008, constituting the first clear example of a 
class 0 YSO (a protostar) with a large accretion burst. The usual assumption that in YSOs accretion 
and ejection follow each other in time needs to be tested. Radio jets at centimeter wavelengths are 
often the only way of tracing the jets from embedded protostars. We searched the Very Large Array 
archive for the available observations of the radio counterpart of HOPS 383. The data show that 
the radio flux of HOPS 383 varies only mildly from January 1998 to December 2014, staying at the 
level of ~ 200 to 300 ^Jy in the X band (~ 9 GHz), with a typical uncertainty of 10 to 20 p,Jy in 
each measurement. We interpret the absence of a radio burst as suggesting that accretion and ejection 
enhancements do not follow each other in time, at least not within timescales shorter than a few years. 
Time monitoring of more objects and specific predictions from simulations are needed to clarify the 
details of the connection betwen accretion and jets/winds in YSOs. 

Subject headings: stars: formation - stars: protostars - stars: radio continuum - ISM: jets and 
outflows 


1. INTRODUCTION 

The physical mechanisms of gas accretion and ejection 
are fundamental aspects in the paradigm of low-mass star 
formation. In its most commonly assumed form, this 
paradigm requires t he quasi-stationa ry evolution of the 
forming stars (e.g., Shu et al. 19941. This assumption 
is seriously challenged by the observational evidence of 
episodic accretion. First, a handful of individual bursting 
objects such as FU Ori and EX Lup provide evidence that 
large accretion bursts in Young Stellar Objects (YSOs) 
do actually ha ppen, with durations from a few months to 
100 yr (e.g., Hartmann fc Kenyon||1996 iHerbig 2008 
Audard et al. 2014|). Second, surveys show that YSOs are 
systematically underluminous with res pect to the expec¬ 
tation of steady-accre tion models (e.g.,|Evans et al.|2009 
Kryukova et al.|2012 ). Episodes of signmcantly increasec 
accretion would allevia te this discrepancy (e.g., Vorobyov 
2009| |Zhu et al.||2009|) . Third, episodic accretion m pro- 
tost ars^^'bFcIassirand I YSOs — has also been invoked 
to explain the observed spread of the Herzprung-R ussell 
diagram of young star clusters ( Baraffe et al.||2012 ). 

Although time monitoring campaigns are revealing a 
large diversity of variation phenomena i n YSOs (e.g., 
Costigan et al.||2014 Gunther et al. 2014), it is still not 
clear whether or not accretion bursts are a widespread 
phenomenon in YSO evolution. The evidence is partic¬ 
ularly scarce for the case of protostars. Only a few class 
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1 or young class 11 objects have been r eported to show 
infrared bursts related to accretion (e.g., Reipurth fc As- 


pin|2004 Caratti o Garatti et al.|20TTl Covey et al.|2l)ll[ 


Fischer et al. 2012| IMuzerolle et al. 1120131). 

Recently, ISatron et al.| (|2U15|) reported an impressive 
a obiec 


Recently, 

outburst in the embedded object HOPS 383 in Orion. 
This object was clas sified as a protostar in the Spitzer 
survey presented in Megeath et al. ( 2012) and in th e 
Herschel+APFiX survey presented in Stutz et al. (2013). 
Based on the ratio of submm-to-bolometric luminosity 
and in the charac teristic temperature of the spectral en¬ 
ergy distribution, Safron et al. (2015) conclude that the 
object is a class 0 YSO, making HOPS 383 the youngest 
protostar known with an accretion-related burst. The 
outburst is most notably seen at a wavelength of 24 /im, 
in which the brightness increased by a factor x35 from 
2004 to 2008. From their extensive IR to submm follow 


Safron et al. 


(2015) find no evidence for significant 


up, _ 

fading of HOPS (>83 up through 2012. 

Radio observations are a key tool to study YSOs: 
the thermal free-free emitting radio jets are often the 
only way to probe the bipolar ejection of material close 
to the pro tostar in the optically-obscured class 0 ob¬ 
jects (e.g., Rodriguez 1997 Anglada et al. 1998). Ra¬ 
dio observations are also useful to probe non-thermal 
(gyrojsynchr otron radiation from active YSO magneto¬ 
spheres (e.g.,|Gudelfl2002||Forbrich et al.|2007|. System¬ 
atic studies of YSOs show that emission at ~ 3.5 cm and 
its variability in class 0 and 1 protostars is dominated 
by the free-free radio jets, while in more evolved class 11 
and HI YSOs this emission is dominat ed by the magne- 


tospheric (gyro) synchrotron emission (Dzib et al. 2013 
Liu et al. 2014), with possible contributions f roiii disk 


photoevap oration (Galvan-Madr id et al. 2014) or weak 
radio jets ( Rodriguez et al.||2014 ). 

In the standard model of star formation, the magne- 
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tohydrodynamical (MHD) launching of jets provides for 
part of the necessary release of specific angular momen¬ 
tum. The amount of ejected material is usually take n to 
be a constant fraction of the accre ted material (e.g., Pu- 
dritz et al. ]2007l |Shang et al.|2007[ ). There is some obser- 
vational evidence tor this assumption to hol d in samples 
of rel atively massive, optically visible YSOs (Frank et al. 
2014 and references therein), but it has yet to be tested 
in the most embedded objects and in time domain. A 
good way of doing this test in embedded class 0 objects 
is to look for possible correlations between the free-free 
emission from radio jets and accretion signatures as seen 
in the near to mid infrared. Bursting protostars offer a 
unique opportunity to perform this test under different 
conditions in a given object. 

In this Letter, we report on pre- to post-burst centime¬ 
ter wavelength observations of the bursting class 0 YSO 
HOPS 383. 


2. DATA SETS AND BASIC RESULTS 

To analyze the time behavior of HOPS 383 at centime¬ 
ter wavelengths, we looked for the available observations 
in the Very Large Arraj0 archive. In none of these ob¬ 
servations is HOPS 383 at the phase center, however, 
it falls inside the field of view, allowing the determina¬ 
tion of its parameters after correction for the primary 
beam response. In the following subsections we describe 
the observation epochs. Table 1 summarizes the relevant 
properties of the data. 

In all the epochs the position of the radio source 
matches within 0.3” to 1.7” with the infrared position re¬ 
ported bylSafron et al.| (|2015l): a(J2000) = 5'*35'"29f81, 
^(J2000) = —4""59'51'.'l. I'his is a fraction of the synthe¬ 
sized beam in all the epochs in which the radio source 
appears unresolved (see Table 1), and about the synthe¬ 
sized beam size for the 2011 epoch with subarcsecond 
angular resolutioij^ Since the next closest radio detec¬ 
tion is « 120” away, and the Spiteer-MIPS angular reso¬ 
lution is ^ 6”, we conclude that the radio detection here 
reported i s the counterpart of the bursting protostar re¬ 
ported by ISafron et al. (12015 1. 


2.1. 1998 

These observations were made in 1998 January 13 as 
part of the project AR 387. The analysis ha s been pre¬ 
sented and discussed in|Reipurth et al.|(|1999|). The data 
were obtai ned at a frequen c y of 8 .40 GHz in the D con¬ 
figuration. iReipurth et al.| ( 1999|) do not report a radio 
counterpart to HUPS 383, most probably because they 
adopted a stringent limit of 5-cr to consider a detection. 
Our reanalysis of these data shows a faint source with 
flux density of 0.23 ± 0.05 mJy that coincides spatially 
with HOPS 383. 


2.2. 2008 

These observations were made in 2008 March 14 and 
June 1 in the C and DnC configurations at 4.86 and 8.46 
GHz as part of project AT359. The data from the two 

® The National Radio Astronomy Observatory is operated by 
Associated Universities, Inc. under cooperative agreement with 
the National Science Foundation. 

^ The distance to the IR position in this epoch is 0.5” to the 
centroid of the radio emission, and 0.3” to the brightest SE peak. 


epochs were concatenated for increased sensitivy, but a 
counterpart to HOPS 383 was not detected with a 4-cr 
upper limit of 0.22 and 0.19 mJy at 4.86 and 8.46 GHz, 
respectively. 


2.3. 2011 

These observations were made in 2011 August 13 with 
the upgraded Jansky Very Large Array (JVLA) in the 
A configuration under project llA-220. The tuning is 
centered at 4.81 and 7.42 GHz, with total bandwidths 
of approximately 1 GHz. This is the only epoch with 
subarcsecond angular resolution. The source associated 
with HOPS 383 is for the first time resolved into two 
components (SE and NW) with a separation of « 0"45 
and joined by a ridge of faint emission. In maps done 
with (u, v) weighting close to natural, we measure flux 
densities that are consistent with rest of the detections 
at other epochs: 0.27±0.02 at 4.81 GHz and 0.31±0.02 
mJy at 7.36 GHz (see Table 1). The spectral index be¬ 
tween these freqnencies is 0.3, indicating free-free emis¬ 
sion with low optical depth. Figure 1 shows a map at 7.36 
GHz done with weighting between natural and uniform 
to maximize angular resolution. The double nature of 
the radio source is clearly seen, probably from a bipolar 
radio jet or a binary. This map with the highest angular 
resolution already suffers from significant (^ 50%) flux 
filtering of extended emission. 


2.4. 2014 

The five more recent observations were made with the 
upgraded JVLA centered at a rest frequency of 9.81 GHz 
(3.0 cm) during 2014 October-December. These obser¬ 
vations are part of project 14B-230. At that time the 
array was in its C configuration. The total bandwidth 
of the continuum observations was about 4.0 GHz. The 
individual flux densities determined for the five epochs 
are listed in Table 1, as well as the flux density from 
the concatenated data. There appear to be marginal 
flux variations within 2014. The maximal flux density is 
293 ± 13 ^Jy on December 12 2014, which is 71 ± 18 /rJy 
above the minimal value of 222 ± 12 /rJy on October 27 
2014. Then, without consi dering the few % uncert ainty 
in the absolute flux scale (Perley & Butler 20131, the 
largest flux difference within 2014 is significant at 4 (t. 
We also tested for intra-day variability by splitting each 
of the five observing epochs into two time chunks. In 
four of the five epochs, the flux measured in each halve 
of the data is within la of the flux reported in Table 1. 
In the November 3 2014 epoch, the fluxes in the halved 
data are within 2a of the epoch average, and this is the 
noisiest of the 2014 observations. We discard significant 
intra-day variability. These results are consistent with 
the relatively long variation timescale expected for the 
radio em ission of class 0 YSOs, compared to class H and 
class HI (Liu et al.||20I4). Finally, the high sensitivity of 
these JVLA observations allowed the data to be split into 
smaller frequency chunks. We made images from the con¬ 
catenated data in four windows, each 1-GHz wide. The 
intra-band spectral index of HOPS 383 is consistent with 
0, indicating optically-thin free-free emission. 

3. DISCUSSION 

3.1. No radio counterpart to the infrared hurst. 
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Table 1 

Data Summary 


Epoch 

Frequency 

[GHz] 

Array, HPBW, PA 
[; arcsecxarcsec; deg] 

Phase Center 
[h:m:s; deg:arcmin:arcsec] 

HOPS383 offset 
[arcsec] 

rms noisd^l 
[/ijy beam^] 

Flux Densitjj^ 

[/^Jy] 

15.01.1998 

8.46 

D; 9.2 X 7.5, -17.8 

05:35:24.395; 

-05:01:07.27 

111.1 

42 

230 ± 50 

14.03.2008+01.06.2008 

4.86 

DnC+C; 6.0 x 4.2, 49.8 

05:35:24.400; 

-05:01:43.25 

138.2 

56 

< 224 

14.03.2008+01.06.2008 

8.46 

DnC+C; 4.8 X 3.1, 57.4 

05:35:24.413; 

-05:00:07.25 

82.2 

47 

< 188 

13.08.2011 

4.87 

A; 0.64 X 0.42, -4.6 

05:35:23.420; 

-05:01:30.52 

137.8 

10 

270 ± 20 

13.08.2011 

7.42 

A; 0.44 X 0.27, -3.6 

05:35:23.420; 

-05:01:30.52 

137.8 

11 

310 ± 20 

27.10.2014 

10.0 

C; 2.18 X 1.66, -1.9 

05:35:23.480; 

-05:01:32.30 

138.5 

13 

222 ± 12 

03.11.2014 

10.0 

C; 2.04 X 1.81, -7.2 

05:35:23.480; 

-05:01:32.30 

138.5 

15 

267+ 16 

06.12.2014 

10.0 

C; 2.29 X 1.70, -28.3 

05:35:23.480; 

-05:01:32.30 

138.5 

11 

293 ± 13 

13.12.2014 

10.0 

C; 2.35 X 1.61, -11.1 

05:35:23.480; 

-05:01:32.30 

138.5 

13 

289 ± 15 

15.12.2014 

10.0 

C; 2.28 X 1.65, -29.3 

05:35:23.480; 

-05:01:32.30 

138.5 

10 

246 ± 10 

All 2014 

9.81 

C, 2.32 X 1.62, -14.1 

05:35:23.480; 

-05:01:32.30 

138.5 

6 

260 + 6 


^ rms noise corrected for primary-beam attenuation around the position of HOPS 383: a = 5^35"*29?81, & = 
Flux density of HOPS 383 corrected for primary-beam attenuation. The errors are o nly statistical, resultini 
the quality of the fit. The absolute uncertainty of the (J)VLA flux scale is a few percent (|Perley Butler|2013| 
are 4 (t. 


from the noise in the maps and 
Upper limits in the 2008 epoch 



Figure 1. High-angular resolution JVLA image of the radio 
counterpart of HOPS 383, observed on August 13, 2011, and 
centered at 7.42 GHz (4.0 cm). HPBW= 0”26 x 0!'21, P.A.= 
—9.0°. The target clearly has two components. The peak 
insensities of the SE and NW components are 75 /rJy beam“^ 
and 55 pJy beam“^, respectively. The rms noise is 11 /rJy 
beam“^. 


MHD-launched jets that load with them a fraction of 
the accreted mat erial are a basic feature of star forma - 
tion models (e.g., Pudritz et al.|2007 Shang et al.|2007 ). 
Furthermore, there is increasing evidence that this is the 
case even for the you ngest systems: the class 0 protostars 
(e.g., Li et al.|j2014 ). In these deeply embedded objects, 
free-free emission from radio jets is one of the best ways 
to tr ace the outflowing material in the inn er few hundred 
AU ( Rodriguezjl997 Anglada et al.|199'8 1. The origin of 
the radio emission i s often considered to be shock-ionized 
gas ( Anglada]|1996 ), but models that include X-ray ion- 
ization within the X-wind scen ario also reproduc e the 
observed properties of radio jets (|Shang et al.|2004|). Re¬ 
gardless of the detailed jet emission mechanism, in the 
models it is assumed that an increase of accretion is fol¬ 
lowed by an ejection enhancement. Therefo re, if the mid- 


IR bu rst truly traces a burst of accretion (Safron et al. 


20151, then a significant increase in the radio tlux should 


follow. 

We find that in HOPS 383 the infrared burst does not 
have a counterpart in the radio. Figure 2 shows the radio 
light curves. The long-term (decadal) radio flux densities 
show mild variations, staying at a level of ^ 200 to 300 
/rJy. The earliest 1998 flux density is the same within 
Icr with the latest 2014 measurements (Table 1). The 
average flux density of the detections in the X band (8 
to 10 GHz) is Sx = 265 /rJy, with a dispersion of 28 
/rJy, or about the noise in single epochs. Taken at face 
value, it may even seem that during the IR-burst the 
radio flux has a minimum (see Fig. 2). We constrain 
the signihcance of a possible radio flux decrease in 2008 
by considering the non-detection in this epoch to have 
a flux between 0 and the 4 (t upper limit. The possible 
flux decrease between the 1998 detection and the 2008 
non-detection has a signihcance between 3.3(T and 0.6(7. 
The respective hux increase between 2008 and 2011 (the 
epoch with the maximum hux) has a signihcance between 
6.0(7 y 2.3(7. These estimates do not include the absolute 
hux uncertainty of a few percent. In the following sub¬ 
sections, we discuss several possible explanations for our 
result. These possibilities do not necessarily exclude each 
other. 


3.2. Protostellar accretion and ejection variations are 
not correlated in general. 

Accretion rates and jet/wind ejection rates have been 
shown to be correlated within samples of the rela tively 
evolved class H YSOs and FU Orionis objects (e.g., Har- 
tigan et al.||1995| |Calvet||1998|). However, the idea that 
in a given object accretion and ejection follow each other 
in time has only started to be tested. Contrary to the 
simplest expectations, it may be the case th at they do 
not follow each other. |Ellerbroek et al. (2014) could not 
establish a relation between outflow and accretion vari¬ 
ability in the Herbig Ae/Be star HD 163296, the former 
being measured from proper motions and radial velocities 
of the jet knots, whereas the latter was measured from 
near-infrared photom e tric an d Bry variability. Similarly, 
Connelley & Greene (2014) monitored a sample of 19 
embedded (class 1) Y SUs with near-IR spectroscopy and 
found that, on average, accretion tracers such as Bry are 
not correlated in time to wind tracers such as the H 2 
and [Fe H] lines. The non-detection of a radio burst in 
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Figure 2. Radio light curve of HOPS 383. The red dots and 
inverted triangles are X-band (8 to 10 GHz) measurements 
and upper limits, respectively. The blue symbols are in C- 
band (4 to 5 GHz). Error bars are Tier. See Table 1 for 
details of the observat ions. The vertical gray stripe marks 
the period of time that Safron et al. (20151 estimated for the 
occurrence of the infrared burst. 

HOPS 383 counterpart to the mid-IR burst is consistent 
with the above mentioned observational results. 

Numerical simulations that include the time-variable 
evolution of protostars can give insight on the possi¬ 
ble time correlation of accretion and ejection . Some 3D 
MHD simulations (e.g., Romanova et al.|2009 ) have been 
performed with enough detail to follow the causal connec¬ 
tion between protostellar accretion and jets-|-winds, but 
they follow the evolution of the system for timescales that 
are somewhat short compared to our observations and to 
FU-Ori like bursts in general. Rather large, long bursts 
such as that observed in HOPS 383 are more likely the re¬ 
sult of such processes as gravitational instab ilities within 
the massiv e young circumstellar disk (e.g., Vorobyov & 
Basu|[2M5 ), but these simulations do not yet resolve the 
physics of jet/wind production. Further work needs to 
be done to simulate what happens to the outflows in pro¬ 
tostars with large accretion bursts. 

3.3. Protostellar accretion and ejection variations are 
not correlated in large bursts. 

Jets are launched through MHD processes that require 
the presence of an ordered mag netic field at the stellar- 
radii scale of the magnetosphere (Shu et al.|1994 | and/or 
further out to a significant fraction of the disk (Konigl & 
Pudritz 20001. If in a bursting YSO the accretion rate be¬ 
comes too high, it is possible that the configuration of the 
magnetosphere and the disk magnetic field become dis¬ 
rupted suc h that the launchi ng of a collimated jet/wind 
is shut off (Hartmann 1998). This ‘jet-quenching’ idea 
needs to be explored in detail in simulations. 

3.4. The interpretation of the tracers. 

The final possibility that we consider is that either or 
both of the assumptions that the radio continuum traces 
free-free emission from a jet and that the mid-IR burst 
is due to an enhancement of accretion are incorrect. We 
argue that this is not the case. (Gyro)synchrotron emis¬ 


2422 

tion from this source 


sion from magnetospheres is known to contribute signif- 
icantly in the more evolved class H and HI YSOs (e.g., 
Forbrich et al. 2007), but in class I and especially class 
0 YbUs the radio emission is most likely dom inated by 
the jet (e.g., Dzib et ar]|2013 Liu et al. 2014). Further¬ 
more, the spectral indices measured in the most sensitive 
epochs are consistent with free-free emission (see sections 
2.3 and 2.4). Regarding the mid-IR burst, altho ugh the 
larges t flux increase is at a wavelength of 24 /im, [Safron 
et al. (2015) showed the occurrence of brightening from 
the near-lR to the submillimeter, making the case for 
a large (~ x35) increase in bolometric - and therefore 
accretion - luminosity. 

A different instance of a misinterpretation of the obser¬ 
vations would happen if the mid-IR burst and the radio 
emission arise from the different components of a binary 
system. Indeed, from Fig. 1, the radio emission can be 
interpreted as either a bipolar jet or a binary, with the 
feature s separated by » 19 0 AU (0.45” at a distance of 
414 pc, Menten et al.||2007). However, even if the target 
is a binary, the lack of a large increase in the combined 
(unresolved) radio flux at any epoch remains to be ex¬ 
plained. 

Finally, it is also possible that there was an ejection 
event traced by a radio jet (as suggested by the morphol¬ 
ogy in the epoch that resolved the source, see Fig. 1), but 
that the total flux density of the jet features did not have 
large variations. Based on a measurement of the proper 
motio ns in the lobes of the radio jet in IRAS 16293- 
Pech et al. (20101) inferred a recent bipolar ejec- 
~lnterestingly, the total flux den¬ 
sity at 8.5 GHz added over the lobes reported by those 
authors only increased by ^ 40% from 2003 to 2008. A 
large increase in the radio flux could rapidly fade if the 
ionized material recombines without further ionization. 
This could happen if some material is dense en ough and 
gets shadowed from the ionizing sourcef^ (e.g., Galvan- 
Madrid et al. 2011). If the possibility discussed here is 
correct, then time monitoring should be done at high 
angular resolution. 

4. CONCLUSIONS 

Our search in the VLA archive for a cm radio counter¬ 
part to the infrared burst recently reported in the class 
0 YSO HOPS 383 yielded a negative result. The lack of 
a counterpart in the radio jet to the accretion burst sug¬ 
gests that accretion and ejection variations do not follow 
each other in time, at least not for large (~ x35) accre¬ 
tion enhancements and within short (up to a few years) 
periods of time. 

Our observations are consistent with recent reports, 
using different techniques, of a lack of temporal corre¬ 
lation between accretion and jet/w ind tracers in class 1 
YSO s (jGonnelley fc Greene 2014) and a Herbig Ae/Be 
star (Fllerbroek et aLpUldj)! We discussed possible in- 
terpretations to these observations in the context of the 
available models of accretion/ejection from th e boundary 
of a stellar magnetosphere with a disk (e.g., Romanova 
et al.j 2009). Time monitoring of more sources and more 
specinc model predictions are needed to clarify the de¬ 
tails of the connection between accretion and jet/wind 

^ The recombination timescale for gas with electron density rie ~ 
10® cm“® is only one month. 
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ejection in YSOs. 
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